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We discuss a method to extract the CKM 7 angle combining Dalitz plot analysis of B — * 
" and untagged B° , B° -» A s 7r + 7r . The method also allows obtaining the ratio and phase 
difference between the tree and penguin contributions from B° and B° — » if* ± 7r =F decays and direct 
CP asymmetry between _B° and B°. From Monte Carlo studies of 100K events for the neutral 
mesons, we show the possibility of measuring 7 with a precision of ~ 5°. 



I. INTRODUCTION 

In the Standard Model, CP violation is implemented 
by the existence of CKM complex parameters. Interfer- 
ence between processes with different weak phases con- 
tributing to a same final state, can generate asymmetry 
on charge conjugate B meson decays, allowing one to 
measure the weak phase difference since the asymmetry 
effects are in some way proportional to it. The common 
method 0, Q to extract the CKM (3 phase, explores the 
interference between the phases generated in the B° — B° 
oscillation, in the decay to the same J/^K S final state. 
For CKM 7, the most well established methods 0, 0, Q 
explores the interference in the time-independent decays 
B~ -> D°K- and B~ -> D°K~, when D° and D° de- 
cays to the same final state. Using these methods, 7 
should be determined at the LHC [6| with precision of 
~ 5°, for one year of data taking. 

Beyond two body interferences, one could explore the 
asymmetry associated with the Dalitz interference be- 
tween intermediary states in three body B decays. This 
was initially proposed for B ± — > n ± TT + n~ 0], where XcO 
plays a fundamental role as reference channel. However, 
the method is statistically limited by the low contribution 
of the Cabibbo suppressed amplitude B ± — ► XcO' r± , es- 
timated on the already observed Cabibbo allowed B ± — » 
XcoK^ channel HQ. Recently, new effort is being made 
to explore 7 in B — > Knir decays. While some approaches 
requires time-dependent analysis [l(| E3 , we first showed 
[12| | how to proceed with an untagged analysis, with the 
benefit of higher statistics. Computer simulations (l3j 
showed that tagging introduces inefficiencies in the order 
of 90% in LHCb data. 



II. METHOD 



The major intermediate sates for B^- — > K ± tt + tt and 
B°, B° -> ifs7r+7r-@,H[Il], are resumed in tableUwith 
the respective tree and penguin contributions, followed 
by CKM weak phases contained in the amplitude. Since 
magnitudes and phases extracted by the fit procedure are 



the overall contributions, it's not possible to isolate and 
measure 7 in a simple amplitude analysis of B + and B~ 
and B° and B°. 





resonance 


contribution 


weak phase 


B+ 


I-r(890) u 7r+ 






A"*(1430) u 7r+ 


K + P {770) u 


V b tVt 3 P + V bu V* s T§ 


7 


A'+/„(980) 




v bc v; s t s 




B° 


K*{890) + tt- 


v bt v t * s p + v bu v; s T c 


7 


A*(1430) + 7r~ 


K S p(770)" 


v H v: s p+v bu v: s T S 


7 


K S f (980) 


K s XcO 


v bc v; s t s 





TABLE I: Resonances and dominant contributions from B — > 
Ktt-k decays . We denote allowed, suppressed by color and 
penguin amplitudes by T c , Ts and P. 



Based in SU(3) flavor symmetry, we expect the same 
penguin amplitudes for the four processes B — » K*tt^ 



and B°,B° — > K* ± it T . The method consist to extract 
the B ± -> K* ir^ penguin parameters 1 , introduce them 
in the B° and B° amplitudes, then use a new technique 
of joint fit to extract the tree phases, where 7 can be 
obtained from. 

As the magnitudes and phases are measured relative 
to a fixed resonance, the compatibility between penguin 
parameters from B ± and B°,B° to K*n, requires that 
the amplitude analysis be made relative to an anchor res- 
onance that have the same amplitude for B charged and 
neutral. We use the non CP violating B — > K XcO am- 
plitude. The asymmetry measured by Belle [l5| for the 
channel B ± AT^ci is A CP = -0.01 ± 0.03 ± 0.02, 
indicating as expected, that the dominant contribution 



1 For simplicity we refer only to K* , while in practice the method 
uses simultaneously the parameters of the resonances if* (890) 
and A"* (1430). 
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is a tree diagram without weak phase. The equality be- 
tween B charged and neutral amplitudes, is based in fla- 
vor SU(3) symmetry considerations. 

Following, there is a schematic representation of the 
method, where arrows represent quantities extracted by 
fit and 6 ± = (5 T ±j): 



B ± -> JTV* 

B° -» K* + tt- 
B° -» K*-tt+ 



: ape 



i §t 



a P e iSp +a T e l9+ 



ap e + ay e 
7 = 7T~ 



(1) 



The method is based in three basic and well accepted 
hypothesis that can be tested during the analysis. First: 
the dominant contribution of -B 
This has been partially confirmed by BaBar 
test, is to check if the B + — » K*ir 
amplitudes are the same. Second 



K*^ is V bt V t * P. 



The 

f and B~ K*tt~ 
the penguin com- 



arc 



ponents from B ± -> K*^ and B°,B° -> K 
equal [H, [l7|- Third: XcO have the same amplitude for 
B ± ^K ± XcO and B°,B° -> K 8 XcO- The experimental 
test for the second and third hypothesis, consists in the 
equality of the tree magnitudes extracted from the inter- 
mediary process B° and B° to K* ± n T . The confidence 
in the result of 7 is based in the fulfillment of all hy- 
pothesis, whereas the failure of one implies in interesting 
unexpected effects. 



III. AMPLITUDE ANALYSIS 

To extract the B ± — » A'+tt+tt - parameters, we apply 
a maximum likelihood fit to the isobaric amplitudes 



A d 



(2) 



where .Ai are functions of Dalitz variables modeled by 
Breit-Wigner distributions times angular functions and 
form factors; i = A"*(890)°, K*(U30)°, p(770)°, / (980) 
and a x , 5 X parameters are fixed. 

Regarding the neutral system, the separation between 
B° and B a is not obvious, having in sight the final state 
Ksk + it~ . In addition, there is mixing among B° and B°, 
introducing time dependence in the probabilities. One 
possibility is to use tagging for creating two separated 
sets of events and apply a maximum likelihood fit using 

M(At) = e -(r/2-*M)At ^ C0S (AmAi/2) 

-iq/pAsm(AmAt/2)], (3) 

M(At) = e -(r/2-iM)A i ^ cos (Am At/2) 

-ip/qA sin (Am At/2)}, (4) 



where A e A are time-independent amplitudes analogue 
to ©, for B° -» As7r + 7r _ and S° -> K s tt + tt~ decays. 

The probability distribution for the final state 
Ks^ + i^~ , independent of it's origin, is given by the sum 
\M(At)\ 2 + \M(At)\ 2 . It was found [ij, [HI that in the 
case \p/q\ = 1, that this sum displays the interesting 
property of canceling mixing and time dependence terms, 

\M{At)\ 2 + \M(At)\ 2 = e- rt {\A\ 2 + \A\ 2 ). (5) 

For the neutral system, we can use the total normalized 
amplitude 



A + £< anz i5i A | 2 + I dxe iS * A x + £ s &iZ iSi A? 



iVi + N 2 



Ni = I \a x e iSx A x +^2 a i ei5 ' A| 2 ds v ds jk, 
J i 

N 2 = \a x e lS * A X +Y1 ^ elS ' A ^ ds *i ds i k ' ( 6 ) 



As a fundamental step in 7 extraction, we developed 
a novel method of joint fit for extracting the parame- 
ters cii, Si and a~i, Si from the same set of untagged time- 
independent events fl2l |. In this method, we apply a max- 
imum likelihood fit using as PDF the eq.©, were three 
parameters a x , S x and S x are fixed. Despite the fact that 
XcO resonance parameters are equal, the term d x is kept 
free instead of being fixed to the same value of a x . This 
should be interpreted as if we multiplied the amplitude A 
by a scaling parameter s. In this case, a~ x = s a x whereas 
the other parameters a~i contains s also. That allows one 
to measure difference in number of events from B° and 
B°, investigating direct CP asymmetry. In one single 
procedure it's possible to extract in an independent way, 
the parameters from the amplitudes A and A and the 
ratio in number of events, calculated by the ratio of the 
normalizations Ni/N^. 

In the convention (Ks, 7r + , 7r _ ) — > (1,2, 3) for the final 
state particle numbering, the charge conjugation opera- 
tion switches the Dalitz variables S12 «-> S13. The reso- 
nances K* bands are centered in different axis for A and 
A, establishing some sort of signature for the event origin 
B° or B°, as can be seen in figure Q] However, in the un- 
tagged analysis, data supplies a joint Dalitz plot, where 
a generic point have an unknown origin. Although we 
cannot distinguish events, the joint fit can identify two 
different overlapping surfaces of B° and B° due to the 
non-overlapping of at least one interference region from 
both amplitudes. The non-overlapping assures that the 
amplitudes can be explored in an independent way by 
the fitting procedure, guaranteeing the unicity in the fit 
result. 

In the Belle B°,B° -> K s t: + tt~ analysis [3, the au- 
thors fit a untagged sample of Kstttt, using as total am- 
plitude the sum |A(cc,y)| 2 + \A(y, x)\ 2 . The amplitudes 
have same parameters of magnitude and phase, assuming 



FIG. 1: Schematic representation for the stand-alone DP of 
a)B° -> A S 7r + 7r~, b)B° -» K S it + TT~ and in c) the joint DP 
obtained from the untagged analysis. Bands are pictorially 
representing resonances regions, where for simplicity we ig- 
nore angular effects in the distribution. Circles are displaying 
regions of interference. Is the non-overlapping of circles that 
assures the separation between both amplitudes in the joint 
fit analysis. 
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no CP violation and differing only by a Dalitz variable 
exchange. Using the joint fit technique, it's possible to 
extract different parameters to both decays and measure 
experimentally CP violation, giving one step ahead in the 
technical difficulty existent so far. 



IV. FEASIBILITY STUDY 

To investigate the feasibility and the error dimension 
in the joint analysis, we generated_and fitted 100 Monte 
Carlo experiments of 100K B°,B a — > J^s7r + 7r _ events 
(the expected number for one year of data taking in 
LHCb [20]). All resonances are described by Breit- 
Wigner distributions and the total amplitude is given by 
(|6|) . The parameters used in thegeneration were inspired 
in the observed parameters by [3, Q • The Dalitz plot dis- 
tribution for one generated experiment of 100K events is 
displayed in figure^ where the small contribution of the 
resonance \c0 can be seen. 

In table HU column input, we show the generated mag- 
nitudes and phases to each resonance of B° and B°. The 
extracted parameters by the joint fit are shown in the last 
column, where the displayed values are the central values 
of the gaussian distributions of the 100 experiments and 
the errors are the width of the same gaussian distribu- 
tion. One can note that the extracted quantities are in 
agreement with the generated ones and with small errors, 
stating the feasibility of the joint fit method. 

One important issue concerning the 7 extraction, is 
the size of the ratio r = ax/ap and the phase difference 
9 = 5? — Sp from the tree and penguin amplitudes of the 
K* resonance. The theoretical knowledge of these quan- 



FIG. 2: B°, B° -> K 3 n + -K~ events distribution from fast MC 
generated with table ITT1 parameters. 



decay B"/B u 


input 


fit 100K events 


XcK B a /a 
So /So 


0.30/0.30 
3.78/3.78 


fixed/(0.30 ± 0.03) 
fixed/fixed 


A*(890)tt ox/5i 
Si /Si 


1.17/1.30 
0.40/5.98 


1.17 ±0.06 /(1.30 ±0.01) 
(0.41 ±0.08)/(5.99 ± 0.07) 


A"o(1430)tt a 2 /a 2 
S2/S2 


2.45/2.72 
0.375/6.00 


(2.45 ±0.11)/(2.72± 0.13) 
(0.39 ±0.08)/(6.00 ± 0.06) 


p u K s a-i/a-i 

s 3 /s 3 


0.60/0.60 
1.20/1.20 


(0.60 ± 0.04)/(0.60 ± 0.04) 
(1.22 ±0.09)/(1.20 ± 0.07) 


f K s a 4 /a 4 
84/84 


1.03/1.03 
2.30/2.30 


(1.02 ±0.06)/(1.04 ± 0.05) 
(2.30 ± 0.07)/(2.30 ± 0.08) 


\N°{B° -> KsTt + tt-)/N°(B -» K s Tv + ir~) = 0.84 ± 0.12 



TABLE II: Monte Carlo simulation for B°, B° -> K s it + tv~ 
decay. We generated sample with the the parameters at and 
Si for B° and at and Si for B°. The third column display the 
fit results with one hundred samples for each 100K events. 
The last line is the ratio between B°,B° — > K s ir + iv~ number 
of events. 



tities is model dependent. Some groups using factoriza- 
tion approach [21(, arrive to large r and small 9. On the 
other hand, non-factorisable approach for pseudoscalar- 
pseudoscalar B decay [HI, presents an opposite scenario, 
with small r and large 9. The joint fit applied to the real 
data, will be able to define which theoretical approach is 
more adequate, since it is possible to measure the param- 
eters under discussion. To our study, we choose r = 0.45 
and generate the experiments with 7 = 69°. 
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If in the unconstrained analysis, the three hypothesis 
prove to be correct, we can assume that the K* reso- 
nances parameters from B° and B° are given by the fol- 
lowing tree and penguin sum: 

ae lS = a P e lSp +a T e i6+ , 

ae a = a P e l5p +a T e l9 \ (7) 

In this case, we can use the scheme (QJ to extract 7 using 
a constrained fit, where the penguin amplitude is fixed 
and the tree parameters ax and 9 ± are directly extracted 
from K* . Using the parameters from tableQll we measure 
7 = 69° ±5°. 

V. CONCLUSION 

We presented a method to extract the CKM 7 angle us- 
ing a combined Dalitz plot analysis from B ± — * K ± tt + tt^ 
and B°,B° — > K s tt + it~. This approach use three basic 
hypothesis that can be tested before one proceeds to the 
constrained fit. For measuring the B neutral parameters, 
we use a new technique of joint fit that allow us to extract 



in an independent way the amplitudes from two summed 
surfaces in a joint sample of B° + B° untagged events. 
We carried a fast MC study to estimate the errors asso- 
ciated to the joint fit technique. Assuming a relatively 
low statistics for the anchor resonance Xco we obtained 7 
with a 5° error. 

During the analysis we can measure CP violation by 
counting the number of events from B° and B°, informa- 
tion extracted in the joint fit, or exploring Dalitz symme- 
tries as discussed in [l8|. We can also measure the ratio 
and phase difference from tree and penguin amplitudes of 
the K* resonance, defining which theoretical approach, 
factorisable or non-factorisable is more adequate. 

The method presented here is competitive with the 
other approaches to determine the CKM 7 angle [|| and 
needs the high statistics expected for the LHCb experi- 
ment, due to the small contribution of the reference chan- 
nel B — * x°K- However, in the case that the fp (980) res- 
onance is dominated by the ss component [231 ] . or even 
if the ratio between the tree and penguin is negligible 
[HJ|, the B — > / (980)X amplitude could take place of 
the charmonium as a reference channel in the analysis. 
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